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A simplified version of White's Density Matrix Renormalization Group (DMRG) 
algorithm has been used to find the ground state of the free particle on a tight- 
binding lattice. We generalize this algorithm to treat the tight-binding particle in 
an arbitrary potential and to find excited states. We thereby solve a discretized 
version of the single-particle Schrodinger equation, which we can then take to the 
continuum limit. This allows us to obtain very accurate results for the lowest energy 
levels of the quantum harmonic oscillator, anharmonic oscillator and double-well 
potential. We compare the DMRG results thus obtained with those achieved by 
other methods. 

PACS numbers: 03.65.-w, 03.65.Ge 



I. INTRODUCTION 

The Density Matrix Renormalization Group (DMRG) method originated from the 
study of a very simple problem: the quantum behaviour of a single particle on a lattice [0. 
The standard renormalization group (RG) method applied to this problem fails completely 
and the understanding of this failure was the key that led to the DMRG. The DMRG has 
since been sucessfully applied to interacting many-body problems in Condensed Matter 
and Statistical Mechanics and it is very promising in Atomic, Molecular, Nuclear Physics 
and Field Theory (see Refs. for an overall account of the DMRG). In view of all these 
developments, it is interesting to come back to the origin of the DMRG and see how it works 



for a single particle under the action of a potential. This problem is not purely academic 
since it provides a testbed for new ideas and techniques relating to the DMRG in simple 
models provided by Quantum Mechanics (QM). In addition, these simple models can contain 
interesting physics. For example, the instantons of gauge theories have a very nice analog 
in the double-well potential. We also hope that the DMRG will provide an accurate new 
method to solve the Schrodinger equation numerically. 



II. THE PROBLEM 

The purpose of this paper is to present a DMRG algorithm for finding the ground state 
and excited states of a particle confined to the real line — oo < x < oo, and whose dynamics 
is governed by the Hamiltonian H = p 2 + V(x), where p 2 is the kinetic energy and V(x) is the 
potential energy. The first step before we can apply the DMRG to QM is to discretize the 
Schrodinger Hamiltonian H. This can be done by constraining the position of the particle 
to take on the discrete values x n = h(n — ^p-), (n = 1, 2, ... , N), where N is the number 
of allowed sites and h is the lattice spacing. The particle is thus confined to a box of size 
R = xtv — X\ = h(N — 1). After this discretization, the Hamiltonian H = p 2 + V(x) becomes 
the N x N matrix 

2/h 2 + V n n = m 

— 1/h 2 \n — m\ = l (1) 
otherwise , 

where V n = V(x n ). Conversely, one can recover the Schrodinger Hamiltonian from Eq. ([I]) 
by taking the continuum limit h — >• 0, N — > oo with R = N x h kept fixed, and then letting 
R — >• oo. In the free case, V(x) = 0, the discrete Hamiltonian, Eq. ([[]), coincides, up to the 
overall factor 1/h 2 , with the Hamiltonian studied by White and Noack in Ref. [EL 
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III. THE DMRG ALGORITHM 



The problem is to diagonalize the N x N matrix, Eq. (|l|), using DMRG methods. The 
basic strategy is to split the box of length N into a left block Bf with £ sites, two middle 
sites •• and a right block B^_ e _ 2 with N — £ — 2 sites, so that the whole system can be 
represented as Bf • •B^_ e _ 2 . If no truncation is performed, Bf contains all the degrees of 
freedom associated with the £ sites on the left hand side, and similarly B§_ £ _ 2 describes all 
the degrees of freedom on the right hand side. 

The DMRG method is based on the idea that the low-energy properties of the system 
can be described by a few degrees of freedom. In particular, if we are only interested in 
the ground state of the system we may simply consider the blocks Bf and B§_ £ _ 2 to be 
described by a single degree of freedom. The effective Hamiltonian of the system then 
becomes a 4 x 4 matrix whose diagonalization is straightforward. Once this is done, the 
next step is to make a partition of the system either as Bf +1 • •B§_ e _ 3 or as B\_ x • •B§_ i _ l , 
i.e. make the left or right side grow by one site, and repeat the diagonalization. One then 
iterates this procedure, moving the position of the partition, £, leftwards and rightwards 
through the lattice. After several such sweeps through the lattice, the DMRG converges to 
a fixed-point solution which reproduces the ground state of the free tight-binding particle 
to high accuracy. 

A. Superblock Hamiltonian 

This method can be generalized to treat the Ne > 1 lowest energy levels, i.e. Ne — 1 
excited states. The basic idea is that the left and right blocks B L and B R must contain 
Ne degrees of freedom, one for each low energy state. The superblock Hamiltonian H$b is 
therefore a (2N E + 2) x (2N E + 2) matrix given by 
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H, 



SB 



H L -v L 
-v[ h CL -1/h 2 



,t 



V 



(2) 



/ 



-1/h 2 h C R 
-v R H R 

where H L and H R are N E x N E matrices, vi and Vr are ^-component column vectors, 
and hcL and hcR are real numbers. The meaning of these quantities is as follows: Hi is 
the Hamiltonian which describes the interactions inside the block B L , —Vi describes the 
interaction between B L and the site next to it in the superb lock B L • »B R and hcL is the 
Hamiltonian on a single site. The quantities Hr, —Vr and hcR for the right half have an 
analogous meaning. The two terms proportional to —1/h 2 in Eq. (0) come directly from the 
off-diagonal terms in Eq. (]!]). 



The Hamiltonian, Eq. (0), describes the superblock B\ • »B 



R 

N-e-2 



and depends on the 



value of £, i.e. Hsb = Hs B - In particular, the diagonal entries hcL and hcR are given by 



L CL 



h 2 



£+1, 



h 



CR 



h 2 



+ v t 



1+2 



(3) 



Since the blocks B L,R contain A^^ effective sites, Hg B can be defined for I = Ne, Ne + 
l,...,N-N E -2. 



B. DMRG truncation 



The superblock can be used either to enlarge the left block by one site, i.e. B\» — > B'^, 
or the right block. Let us examine how to do this for the left block. For iV^ = 1, the ground 
state of the superblock Hamiltonian, Eq. has four components, which we will designate 
(aL,acL,acR, a R)- The projection of this state onto B L » yields (aL,aci), which must then 
be normalized by dividing by iV a = 
are then given by 



a 2 L + Oql- The new effective Hamiltonian, H' L , and v' T 



(a L , a CL ) 



( 



H L -v L 
-v L h C L 



\ ( „/ \ 



(4) 



y a CL j 



and 



v L = a CL , 



(5) 



where a' L = a L /N a and a' CL = a C L/N a . 

We next describe how to generalize this construction to N E > 1. First we have to obtain 
the lowest Ne eigenstates of H$b, which we shall denote as {(a-L,i, a>cL,i, a>CR,i, a R,i)}f=i, 
where a^i and a^j are A^-component vectors. These vectors are projected onto a set of Ne 
vectors of the block B L », i.e., {(a.L,i,acL,i)}^ v Since this set of vectors is not orthogonal 
in general, we ort honor malize them explicitly. Let us call the set of vectors so obtained 
{(a'i,^, a' CL i )} _ \ The projection B L » — > B' L is performed by an (N E + 1) x N E matrix 



A' 



• a L)Ne 
y a'CL,i ■ ■ ■ a'cL,N E J 



(6) 



Operators associated with B' L are transformed into the new basis via 



H' L = A't 



H L -v L 



-v L h 



CL 



A' 



(7) 



and 



V L,i — a CL,ii 



1,...,^ 



E 



(8) 



Observe that we do not need to construct a density matrix to define a unique projection. 
This is a peculiarity of the single-particle nature of the problem. In a many-body problem, 
we would have to perform the projection using the reduced density matrix due to the non- 
single-valuedness of the projected wave functions. 



C. Initialization and sweeps 

The finite-system DMRG algorithm consists of a warm-up phase in which the system is 
built up to its actual length N, followed by several leftwards and rightwards "finite-system" 
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sweeps which are repeated until convergence is achieved. The initial superblock of the warm- 
up phase is taken to be • »-B^ , where H L (H R ) are taken to be the first (last) N E 
entries of the matrix, Eq. ([I]), i.e. 

H L = (h njm ,n,m = 1, . . . , N E ) (9) 
H R = (hn, m , n, m = N - N E + 1, . . . , N) . 

Similarly, v E = (0, 0, ... , 1) and vr = (1, 0, . . . , 0). For a symmetric potential, i.e. V(—x) = 
V(x), the quantities H r ,vr for the right side can be obtained from the left ones by the 
reflection operation, n — > N E — n + 1. This property applies at each step of the warm- 
up phase of the algorithm and corresponds to the well-known reflection operation of the 
infinite-system DMRG algorithm for interacting systems. Of course, for a non-symmetric 
potential the right block, B R , cannot be formed by reflecting the left one, B L . 

The behavior of the location of the added site, £, during the warm-up process and the 
first sweep can be summarized in the following scheme: 

warm — up : i = N E , . . . , y — 1 
left -> right : i = f , . . . , N - N E - 2 

2 , (io) 

right -> left : i = N — N E — 2, ... , N E 
left -> right : t = N E ,...,%-1 

where we take N to be even. We define a sweep to be a complete cycle that ends when the 
two sites of the superblock B\ • *B^_ t _ 2 are in the middle of the chain. 

IV. ANALYSIS OF ERRORS 

In this paper we shall present the results obtained with the DMRG algorithm applied 
to three well-known potentials and compare them with the results of other methods. The 
three potentials we shall consider are the harmonic oscillator, the anharmonic oscillator and 
the double-well potential. 
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In order to take the continuum limit, we enlarge the number of lattice sites, N, while 
taking the discretization step, h, to be smaller and smaller so that R = N x h. During 
this process, it it important to determine the accuracy of the DMRG results for the energy 
levels as N increases. To this end, we first compare the DMRG results with the known exact 
results for the tight-binding particle with fixed boundary conditions at the ends. The exact 
spectrum is given by: 

%,(:) = n^±±^ 3 , £ <»> = W(||±i>) (ii) 

where n = 0, 1, . . . ,N — 1 denotes the energy levels, j — 1, 2, . . . , N are the lattice sites and 
the N n are normalization constants. 



1 2 3 4 5 
i 1 1 1 1 1 1 1 r 




FIG. 1. Logarithm of the relative error in the ground state energy of a free particle on a 
tight-binding lattice as a function of the logarithm of the number of sites of the chain, varying the 
number of targeted states, Ne- 

We have made runs targeting the Ne = 1,2,3 lowest-lying states and have computed 
the relative error of the energy levels, 

SEi = \{Ei{DMRG) - Et^/E^l i = l,...,N E . (12) 
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For Ne = 1, a slowing down of the convergence and, eventually, numerical instabilities due 
to roundoff error appear on chains of length of order 500 or greater due to a vanishing matrix 
element connecting the blocks. Since this problem does not occur when Ne > 2, we will 
take Ne > 2 in the following. 

In Fig. |TJ, we plot the relative error in the ground-state energy for Ne = 2 and 3 as 
a function of N. These results were obtained with the finite-system algorithm described 
above, taking a sufficient number of finite-system sweeps to obtain convergence. Typically, 
convergence was achieved after five sweeps. In this plot we see that for a fixed number of 
targeted states, Ne, the error increases with the number of sites N. The range goes from 
5E = 5.1 x 10" 15 for iV = 10 to 5E = 1.4 x lO" 10 for N = 60,000 (N E = 2), while 
5E = 2.6 x 10~ 15 for iV = 10 and 6E = 4.4 x 10" 9 for N = 60,000 (N E = 3). Likewise, 
when the number N E is increased we observe that the errors decrease and this effect is more 
apparent when the number of sites N is bigger. In Fig. 0, we plot analogous results for the 
first excited state, and find that the error behaves essentially the same as for the ground 
state. 



1 2 3 4 5 




"i 1 1 1 1 1 1 1 r 

1 2 3 4 5 

Log(N) 



FIG. 2. Logarithm of the relative error in the first excited state energy of a free particle on a 
tight-binding chain as a function of the logarithm of the number of sites of the chain, varying the 
number of targeted states, N E . 
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The lesson to be learned from this analysis is that care must be taken when going to 
large lattice sizes because the DMRG tends to loose numerical accuracy with increasing N. 
On the other hand, targetting a larger number of states N e actually improves the accuracy 
of each state. In going to the continuum limit, there is therefore a tradeoff between the 
discretization error which can be reduced by taking iV large, and the error in the DMRG 
results ||. 

V. THE HARMONIC OSCILLATOR 

We first apply the continuum limit of the DMRG to a simple and exactly solvable model: 
the harmonic oscillator. We write the Hamiltonian 

H = P 2 + X 2 1 (13) 

where P = —i4z- The corresponding exact spectrum with this normalization is E n = 2n+l, 
n — 0, 1, . . . We can use this example to calibrate the accuracy of the DMRG by comparing 
with the exact solution. In Table [I], we present the results for the two lowest eigenstates, for 
which E = 1 and E\ = 3. These results are obtained using five finite-system sweeps. We 
have checked that we obtain the same results to within roundoff error when performing up to 
twenty sweeps. We measure the DMRG energies during the final sweep in the configuration 
in which the left block and the right block are the same size, i.e. the single sites in the 
superblock are at the middle of the chain. The error bars in the table are given by the 
amount of variation in the energy during the last finite-system sweep, i.e. the value given is 
obtained in the last diagonalization step, but the digits in parentheses vary during previous 
diagonalizations in the sweep. We call this method of estimating the error the Global 
Convergence Criterium (GCC). However, the GCC is overly restrictive, as can be seen by 
examining the difference between the continuum ground state energy and the DMRG ground 
state energy, plotted in Fig. | as a function of the DMRG step during five finite-system 
sweeps for N = 3000 and h = 0.01. One can clearly see regions in which the solution has 
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been stabilized as well as depleted regions near the ends of the chains in which the DMRG 
energy is higher. In the inset, we plot the error in energy for a single back-and-forth sweep 
on an expanded scale. This suggests that the appropriate region to measure the energy must 
be located away from the ends. We refer to this as a Local Convergence Criterium (LCC). If 
we adopt the LCC, then the number of stabilized digits increases. This is something which 
we will examine when comparing DMRG results and exact results in the forthcoming tables. 
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FIG. 3. The difference between the ground state energy of the harmonic oscillator obtained with 
the DMRG during 5 sweeps for N=3000 and h=0.01, and the exact continuum value. The data 
plotted starts after convergence in finite-system sweeps has been achieved. 

In Fig. ^ the DMRG ground state energy converges variationally from above to the 
exact ground state energy (note that .Eo(DMRG) is negated in the plot) for a particular 
discretization of the system. This energy is close to the plateaus in Fig. [3] and is lower than 
the exact continuum energy, i.e. there is a discretization (or finite-size) correction which is 
not necessarily variational, and which is negative here. This can also be seen in Table | for 
the exact diagonalization as well as for the DMRG results. 
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Method 


N 


h 


R = N x h 


1.0 -E (N) 


3.0- (AT) 


DMRG 


100 


0.1 


10 


6.253912540(370) x 10~ 4 


3.128521663(378) x 10~ 3 


Exact Diag. 


100 


0.1 


10 


6.253912540 310 x 10~ 4 


3.128521663 376 x 10~ 3 


DMRG 


1,000 


0.01 


10 


6.24989(322) x 10~ 6 


3.1243338(564) x 10~ 5 


Exact Diag. 


1,000 


0.01 


10 


6.24989 256 x 10~ 6 


3.1243337 792 x 10~ 5 


DMRG 


5,000 


0.002 


10 


2.498(36) x 10~ 7 


1.2427(28) x 10~ 6 


DMRG 


10,000 


0.001 


10 


6.2(45) x 10 -8 


3.051(90) x 10~ 7 


DMRG 


20,000 


0.0005 


10 


1.(56) x 10~ 8 


7.0(80) x 10~ 8 


DMRG 


20,000 


0.001 


20 


6.2(64) x 10~ 8 


3.12(50) x 10~ 7 


DMRG 


50,000 


0.0002 


10 


(2.87) x 10~ 9 


(4.99) x 10~ 9 


DMRG 


50,000 


0.0004 


20 


(0.961) x 10~ 8 


(5.00) x 10~ 8 


DMRG 


100,000 


0.0001 


10 


(-0.71) x 10~ 8 


(-4.15) x 10~ 8 


DMRG 


100,000 


0.0002 


20 


(0.16) x 10~ 8 


1(.24) x 10" 8 



TABLE I. The difference between the exact continuum energies, E n = 2n + 1, and the energies 
on the discretized lattice, Eq{N) and E\{N), for the ground and first excited states of the harmonic 
oscillator obtained with the DMRG and exact diagonalization. Here is the number of sites on 
the chain, h is the discretization step and R is the size of the continuum system. The errors (digits 
in parentheses) are determined using the GCC procedure explained in text. 
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In Table | we present several runs in which the final size R of the continuum system is 
varied in order to check that a convergence has been achieved to a certain number of digits. 
If R is large enough, the effect of changing it is not very large due to the exponential fall-off 
of the wave functions. We find similar dependence of the error on R for higher excited states. 
Thus we find that the exact spectrum is very well reproduced by the continuum limit of the 
DMRG. 

In Table | we also show, when possible, energies obtained via exact diagonalization of 
the Hamiltonian matrix, Eq. (|l|). This is possible for system sizes of up to about iV = 2000. 
Exact diagonalization is much more expensive in computer time than DMRG, and yields 
energies for a particular system size that agree with the DMRG to 14 significant digits. 
This level of agreement is quite remarkable given that the largest matrix diagonalized in the 
DMRG procedure is 6 x 6. 

VI. THE ANHARMONIC OSCILLATOR 

Now that we have checked that the continuum limit of the DMRG accurately reproduces 
a well-known solvable case, we apply DMRG to the anharmonic oscillator, which has no 
closed analytical solution (see f| and references therein) . Here our purpose is to compare the 
efficiency of DMRG to that of other standard methods employed in single-particle quantum 
mechanics. We treat the Hamiltonian 

H = P 2 + X 2 + CX 4 , (14) 

where C is a positive coupling constant and we have normalized the mass term. 

In Table [H[ we display the DMRG results obtained using a similar analysis as for the 
harmonic oscillator. We compute the two lowest energy states, compare them with Exact 
Diagonalization results, and obtain the same accuracy as for the harmonic oscillator. We 
also compare with the results of a number of methods commonly used in standard quantum 
mechanics such as the Hill determinant method , Borel-Pade approximants of the pertur- 
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bation series || and variational computations in the energy basis of the harmonic oscillator 

We have computed the state energies for various coupling constants C = 0.1, 1, 10, rang- 
ing from weak- to strong-coupling. We notice that while the Hill determinant method and 
the Borel-Pade perform better for weak couplings, the continuum limit of the DMRG per- 
forms equally well for the whole range of couplings. This is due to the variational nature of 
the method. 



Method 


N 


h 


R = N x h 


E 


Ei 


C 


DMRG 


30,000 


0.001 


30 


1.06528543(310) 


3.30687158(195) 


0.1 


DMRG 


20,000 


0.001 


20 


1.06528543(289) 


3.30687158(195) 


0.1 


Borel-Pade § 








1.06528550 954 




0.1 


Hill meth. § 








1.06528550 954 




0.1 


DMRG 


30,000 


0.001 


30 


1.39235148(046) 


4.64881163(809) 


1 


DMRG 


20,000 


0.001 


20 


1.39235148(038) 


4.64881163(809) 


1 


Borel-Pade § 








1.392350(6) 




1 


Hill meth. § 








1.3925316 




1 


DMRG 


30,000 


0.001 


30 


2.449173(484) 


8.5989993(093) 


10 


DMRG 


20,000 


0.001 


20 


2.449173(484) 


8.5989993(093) 


10 


Variational || 








2.449174 072 




10 


Borel-Pade |] 








2.440(527) 




10 


Hill meth. § 








2.4491740 
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TABLE II. DMRG results for the Anharmonic Oscillator H = P 2 + X 2 + CX A in the continuum 
limit. The notation is the same as that in Table |. 
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From Table J|, one can see that the agreement of the DMRG with the other methods is 
excellent. Moreover, once the DMRG results have converged, the results are lower in energy 
than those of the other methods. Thus, since DMRG is variational, we conclude that the 
DMRG results are closer to the exact results than those of the other methods. 



VII. THE DOUBLE- WELL POTENTIAL 

We hope that it is now clear that the DMRG method is an excellent method to compute 
the energy spectrum and wave functions in quantum mechanical problems. We will now 
apply it to an anharmonic oscillator with a potential in the shape of a double well. This 
problem is particularly interesting for several reasons. One, as we shall see, is that the system 
has a tunable gap which can be used to investigate the dependence of the convergence of the 
DMRG on the size of the gap, and the other is its potential as a new non-perturbative method 
for Quantum Field Theory problems. We shall not pursue the latter goal here but will only 
point out that a number of non-perturbative Renormalization Group techniques have been 



established since Wilson's original formulation of the RG [Kj. These formulations are called 
Exact Renormalization Group (ERG) because they are based in exact RG flow equations 
(see for a recent review on this subject). However, since these exact equations are not 



exactly solvable in general, one usually has to resort to approximate methods in practical 
applications. The question then arises as to how good these approximations are. In order 
to check their validity, they are applied to well-known problems in single-particle Quantum 
Mechanics. The rationale is that if they are not even able to quantitatively reproduce the 
physics of these simple systems, their application to truly field theoretical problems would 
be even less successful. 

The Hamiltonian for this potential reads 

H = P 2 -X 2 + CX 4 , (15) 

where C is the coupling constant and we have normalized the negative mass term. The 
potential has two minima at the positions 

14 



Classically, these minima are degenerate in energy. If treated perturbatively, quantum fluc- 
tuations can modify the classical energy but not lift the degeneracy Splitting of the energy 
levels can only occur if quantum tunneling between the two wells is taken into account. 

In this symmetric potential, the energies can be arranged into pairs E^> depending on 
their even (+) or odd (— ) parity. The energy gap from the ground state to the first excited 
state is then defined as 

A (C)=Et\C)-E { +) (C). (17) 
In the weak coupling limit (C very small), the gap A (C) can be computed using an 



instanton approximation plus higher corrections fll2], giving the asymptotic formula 



A (r) _JV2 g 71 2V2 6299 2V2 2 2 2691107 2^2 

(18) 

This exponentially decreasing behavior produces an essential singularity in the energy gap. 
Our purpose is to capture this highly non-trivial behavior with the DMRG. 

One can understand what happens to the system in the weak coupling limit on physical 
grounds. From Eq. ([16]), we see that the distance between the minima diverges as C — > 0, 
so that we effectively end up with a system formed by two independent potential wells. 
The system then becomes exponentially degenerate for each pair of energy levels and thus 
gapless. 

It is also interesting to use this quantum mechanical example to test the behavior of 
DMRG for gapless systems. This is a very important issue when dealing with the strongly 
interacting quantum many-body systems to which DMRG is usually applied. Since its 
early development, it has been known that the the DMRG method ||T| produces much more 
accurate results for finite correlated systems with a gap than for gapless systems. For gapless 
systems one has to use the finite-system algorithm on larger system sizes than for gapful 
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systems in order to obtain results of comparable accuracy [ 13| . We shall show below that 
the DMRG handles well the case in which the gap between the ground state and the first 
excited state becomes negligibly small, i.e. when the two minima of double well potential 
are far apart. 



Method 


C 


N 


h 


Eq 


A 


DMRG 


1 


1000 


0.01 


0.65764425361(29) 


2.176825710(298) 


Exact Diag. 


1 


1000 


0.01 


0.6576442536117 


2.176825710302 


DMRG 


1 


20,000 


0.0005 


0.657652983(568) 


2.17688305(297) 


R-R 


1 






0.657653005181 


2.17688319705 


DMRG 


0.6 


20,000 


0.0005 


0.39195261(873) 


1.6332847(928) 


R-R 


0.6 






0.39195263337 


1.6332848846 


DMRG 


0.1 


40,000 


0.001 


-1.26549292(138) 


0.11243368(739) 


R-R 


0.1 






-1.26549283721 


0.11243370614 


Inst. 


0.1 








0.11447450849 


DMRG 


0.06 


60,000 


0.0005 


-2.82363949(203) 


0.0072997661(673) 


R-R 


0.06 






-2.82363945845 


0.0072997526870 


Inst. 


0.06 








0.0073139070463 


DMRG 


0.02 


90,000 


0.0004 


-11.106472434(074) 


2.1(074) x 10~ 9 


R-R 


0.02 






-11.106472414954 


2.1043 x 10~ 9 


Inst. 


0.02 








2.10737 x 10" 9 



TABLE III. DMRG results for the double-well potential H = P 2 - X 2 + CX A in the continuum 
limit. The notation is the same as that in Table |. R-R stands for the results obtained with the 
Rayleigh-Ritz method explained in the text. 
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In Table |T|, we present the results for the lowest two energy levels and the gap upon 
varying the coupling constant from C — 1.0 to C — 0.1. These values range from the strong- 
to the weak-coupling regime. The DMRG results are obtained as described previously after 
five finite-system sweeps. We also compare with Exact Diagonalization methods when 
N < 2000 and find excellent agreement, typically up to 10 digits or better. 

In Table |T|, we have also computed the energy variationally using the Rayleigh-Ritz 
method |§ (R-R). We compute the expectation value of the Hamiltonian, Eq. (|l5l), in 
the energy basis of the harmonic oscillator consisting of up to M = 1000 states. In this 
representation, the non-vanishing elements of the Hamiltonian lie within a band and are 
given by 

111 

(n\H\n) = C[-n(n - 1) + ~(2n + l) 2 + -{n + l)(n + 2)] 

Q 

(n\H\n + 2) = -yj(n + 2){n + 1) + C[y/(n + l)(ra + 2)(n + -)] 

(n\H\n + 4) = C-y/(n + 4)(n + 3)(n + 2)(n + 1) . (19) 

In this representation, the Hamiltonian is already in the continuum limit. The Rayleigh- 
Ritz theorem states that the M resulting energy levels will be upper bounds to the first M 
exact energy levels. As seen in Fig. |], the agreement between exact and DMRG methods is 
excellent for this range of the coupling constant and, in fact, the curves appear overlapped 
in the plot. 

Recently, one ERG method has been applied to the study of the double-well potential 
with the aim of probing the method in the whole range of coupling constants f[4| . This 
method is based on the solution of the local potential approximation of the Wegner- Houghton 
equation [pj|,|]I|. The outcome of these investigations is that the ERG performs very well 
in the strong coupling regime where no perturbation treatment is available. However, in 
the weak coupling limit, the ERG fails to reproduce the behavior found with the instanton 
formula, Eq. dl8|). 
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FIG. 4. Energy gap Aq(C) in the anharmonic double-well potential as a function of the coupling 
constant C in the very weak coupling regime (Table HU). Here the results of the three methods: 
the DMRG, Instanton calculation and Rayleigh-Ritz (R-R) overlap. Results for a wider range of 
C (Table |ll|) are shown in the inset. 



To check the performance of the DMRG as compared to the instanton formula, Eq. (|i~8D, 
in the very weak coupling limit, we have extended our computations from C = 0.1 down to 
C = 0.02, where the gap becomes as small O(10~ 9 ). In Table |TT| and in Fig. |], we present 
the results from the DMRG, R-R and the instanton formula. We see that DMRG is an 
excellent non-perturbative method in the entire range of coupling constants and that it is 
able to quantitatively capture the exponentially decreasing behavior of the gap. 



VIII. CONCLUSIONS AND PROSPECTS 

The DMRG method we have presented in this paper is a natural extension of one in- 
troduced by White to illustrate the DMRG algorithm for the most simple problem: a free 
particle on a tight-binding chain ||. We have shown that the DMRG works with very 
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high precision, yielding an accurate determination of the lowest energy levels for three dif- 
ferent potentials: the harmonic oscillator, the anharmonic oscillator and the double-well 
potential. Its performance is better than or comparable to other known perturbative and 
non-perturbative methods. For single-particle Quantum Mechanics, the DMRG does not 
require the use of a density matrix. The number of states retained in the RG procedure 
is equal to the number of states to be obtained, Ne- Aside from Ne, the variational wave 
function of the DMRG has no adjustable parameters. 

Single-particle Quantum Mechanics has been widely used in the past as a testing ground 
for concepts or techniques that can be applied to more complex systems. With this in mind, 
we have studied quantum tunneling through a potential barrier for the double-well potential 
and found a value of the gap very close to the exact one for a large range of coupling 
constants. This is an interesting result because it shows that the DMRG can, in principle, 
handle tunneling phenomena better than other methods such as the Exact Renormalization 
Group [14]]. An interesting topic for future work would be to explore to what extent this 
feature holds for many-body systems or for field theory. 

Although the DMRG was originally developed as a ground state technique, there have 
recently been new developments in using it to obtain dynamical information ]T7],|rS|. In the 
context of single-particle Quantum Mechanics, one could ask whether the DMRG could give 
information about phase shifts, decay rates, etc. These and other questions remain to be 
investigated. 
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